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Work relations for a system governed by Tsallis statistics 


Ian J. Ford and Robert W. Eyre 

Department of Physics and Astronomy, University College London, Gower Street, London WCIE 6BT, U.K. 

We derive analogues of the Jarzynski equality and Crooks relation to characterise the nonequilib¬ 
rium work associated with changes in the spring constant of an overdamped oscillator in a quadrat- 
ically varying spatial temperature profile. The stationary state of such an oscillator is described 
by Tsallis statistics, and the work relations for certain processes may be expressed in terms of q- 
exponentials. We suggest that these identities might be a feature of nonequilibrium processes in 
circumstances where Tsallis distributions are found. 


I. INTRODUCTION 

Work relations are remarkable identities associated 
with the behaviour of systems undergoing nonequilibrium 
thermodynamic processes. As the name suggests, they 
refer to the performance of mechanical work on a system 
brought about by a protocol of change in external force 
fields while the system exchanges heat with a coarsely 
specified environment. The key point is that they hold 
when the thermomechanical processing of a system takes 
place at an arbitrary rate, and not just in the quasistatic 
limit associated with equilibrium thermodynamics. They 
are therefore statements about real nonequilibrium dissi¬ 
pative behaviour and have attracted attention since they 
provide a perspective on the nature of entropy produc¬ 
tion HHH]. 

The identification of principles underlying the transfer 
of energy in the form of heat and work motivated the 
initial development of thermodynamics in the 19th cen¬ 
tury, but the construction of a framework to justify such 
principles from a microscopic perspective is still ongoing. 
Work relations, and the underlying fluctuation relations 
or theorems, rely on rather few assumptions about the 
dynamics of system components and have proved to be 
extremely valuable in this regard. One of the first results 
in this category was developed by Bochkov and Kuzovlev 
[13 E] to describe the statistics of nonconservative work 
performed on a system in an isothermal heat bath. The 
Evans-Searles Fluctuation Theorem [13 US] provided in¬ 
sight into the behaviour of deterministic complex sys¬ 
tems driven by external mechanical and thermal interac¬ 
tions. A wider appreciation of fluctuation theorems, and 
of the significant development that was then underway 
in nonequilibrium statistical physics, followed from the 
derivation of the Jarzynski equality [muH], an elegant 
work relation that can be established from a number of 
perspectives. Further interest in the field was generated 
by the Crooks relation m, a connection between the 
statistics of work performed in two isothermal processes 
that are driven in a mutually time-inverted fashion. Both 
relations require the system to start out in thermal equi¬ 
librium. Experimental relevance of these results has been 
demonstrated |20| . 

We wish to explore work relations that hold for a sys¬ 
tem exposed to an environment with a temperature gra¬ 
dient. The control, characterisation and exploitation of 


nonisothermal conditions at the nanoscale has received 
considerable attention in recent years |21fl23| . Thermo¬ 
dynamic concepts such as system temperature can re¬ 
main valid at these scales [E], though it is clear that 
temporal and spatial fluctuations will be important, as¬ 
pects that can be accommodated through using stochas¬ 
tic thermodynamics as a theoretical description [TU]. The 
basis of this approach is a framework of stochastic system 
dynamics in which the concept of a temperature gradient 
enters through a specification of the statistical properties 
of the environmental noise. 

We have in mind applications to experimental systems 
such as tweezers [23 and nanomechanical devices [26] 
that can be used to manipulate small objects such as 
colloidal particles. There are challenges in establishing 
and controlling a temperature profile at sub-micron spa¬ 
tial scales, but progress has been made in this direction 
|27H29| . Models of heat transport can be used to de¬ 
scribe the elevation of local temperature brought about 
by laser absorption in a liquid medium m and by Joule 
heating in conducting nanostructures m, and it seems 
that a range of thermal profiles can be established using 
suitable arrangements of heat sources and sinks, together 
with choice of media [lain]. 

In such a nonisothermal environment, a system cannot 
achieve thermal equilibrium, and our aim is to under¬ 
stand how the Jarzynski and Crooks results need to be 
adapted. Work relations have been considered for other 
systems that are prevented from reaching thermal equi¬ 
librium, such as glasses [32] • If the dynamics prevent or 
delay evolution towards the true equilibrium, a free en¬ 
ergy of the final state might emerge as a valid concept, 
and so too might work relations. 

We focus our attention on the specific case of a harmon¬ 
ically bound particle in a quadratically varying thermal 
profile. We employ overdamped stochastic dynamics to 
model the behaviour of such a particle showing that is 
governed by Tsallis rather than Boltzmann statistics in 
the stationary state [33 ■ We find that by using a suitable 
notation, the two work relations can be written in a very 
intuitive form for certain classes of protocols. A similar 
generalisation of the work relations to Tsallis statistics 
was proposed by Ponmurugan |34| in the context of a 
basis in Tsallis entropy. In contrast, our results emerge 
from a more standard thermodynamic framework where 
Gibbs-Shannon entropy represents the uncertainty of the 
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state of a stochastic system. We review the nonisother- 
mal stochastic dynamics of systems in Sectionjllj describe 
the Tsallis oscillator in Section |III| and then present ana¬ 
logues of the Jarzynski and Crooks results in Sections 
[TV| and m We consider whether existing experimental 
techniques could be used to establish nonisothermal con¬ 
ditions to a degree that suits our purpose in Section |VI| 
before giving our conclusions in Section |VII[ 

II. DYNAMICAL FRAMEWORK 


possible under the dynamics, and represent either a path 
integration (•) = f dxP\x\{-) where 7^ is a trajectory 
probability density, or the average (•) = / dWp{W){-) 
where p is the probability density for the work done. Eq. 
Q is the celebrated Jarzynski equality, and it holds for 
an arbitrary work protocol as long as the system is ini¬ 
tially in thermal equilibrium. 

Under similar isothermal conditions and initial equi¬ 
librium, the Crooks relation 

Pf{W) =p_R(-W)exp[(W-AF)/A: bTo], (5) 


Dynamical approaches that do not introduce noise are 
available EE], but a stochastic framework jSSHSTj has 
a number of intuitive and mathematical advantages m- 
We consider a system that is weakly coupled to a com¬ 
plex environment where the uncertainty in the interac¬ 
tions that take place between them as time passes is rep¬ 
resented by a stochastic force in the system dynamics. 

We focus attention on a particle that moves in one 
spatial dimension under the influence of a confining po¬ 
tential, a dissipative force proportional to the particle 
velocity and a stochastic force with a strength related 
to the local temperature of the environment but lacking 
autocorrelation. The dynamics in the overdamped limit 
are specified by 


my ox \ my J 

where x is the particle position, m is its mass, (j) is the 
potential, y is the friction coefficient, ks is the Boltz¬ 
mann constant and is the temperature of the environ¬ 
ment. The noise ^ has statistical properties = 0 

and = S{t — t') where the brackets denote an 

average over realisations. Equivalently, we can express 
the dynamics using the following Itb-rules stochastic dif¬ 
ferential equation 


can be derived, where the suffices F and R refer to pro¬ 
cesses that are time-inverted partners, in that the po¬ 
tential evolves over time in process R according to the 
reverse of the evolution in process F. The free energy 
difference in Eq. (§ refers to the outcome of the F pro¬ 
cess. 

Our purpose is to derive analogues of Eqs. Q and 
© for work processes that take place in a nonisother¬ 
mal environment in the limit of overdamped dynamics. 
The appropriate form of stochastic dynamics in these cir¬ 
cumstances has recently received some attention |!MlJ9| . 
particularly with regard to thermodynamic differences 
arising from the neglect of housekeeping entropy produc¬ 
tion in an overdamped treatment gni. The key point 
we exploit is that with overdamped dynamics, the en¬ 
tropy production in the stationary nonisothermal state 
vanishes, such that it takes on something of the charac¬ 
ter of an equilibrium state, and this analogy allows us to 
construct an effective free energy change associated with 
a process and hence to derive work relations. 

The approach we take is to recast the dynamics of Eq. 
Q in an equivalent isothermal form 


1 d^{y,t) 
my dy 


2A:bTo 

my 


1 
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( 6 ) 


dx = - dt + 

my 


/ 2kBTr{x) 
\ my 


dW, 


( 2 ) 


where (f>' = d<j)/dx and dW is an increment in a Wiener 
process with {dW) = 0 and ((dW)^) = dt. 

The work performed on the system as it follows a tra¬ 
jectory (denoted x) in an interval 0 < t < t is 

= £ W{x{t),t)dt, (3) 


in terms of a stochastic variable y{x) and an effective 
potential <i>(j/, t), and to reinterpret the work relations 
that emerge in this representation. The mathematical 
reformulation is straightforward. Using Ito calculus we 
write 


dy 


dy ksTr d^y \ dy 

my dx rri'-f dx^ J dx 


/ 2kBTr 
\ m7 


1 

2 


dW,(7) 


and we choose y{x) such that dy/dx = (To/Tr)^/^ and 


and for a constant environmental temperature = Tq 
it can readily be demonstrated iniini that 

(exp (-lU/fcsdd)) = exp(-AE/A:BTo), (4) 

where AF = F{(jF) — F{(jP) is the difference in Helmholtz 
free energy, at temperature Tq, brought about by a 
change in the potential from (jp = (/)(x, 0) to (jf = (f){x, r). 
The brackets denote an expectation over all trajectories 


d^ 

dx^ 


rpll2rp/ 
1q Ij. 


2t; 


■ 3/2 


( 8 ) 


where T^. = dTj./dx, giving 


\ 7717 zm7 / \lr 


2k^Y 

my ) 


dW, 


( 9 ) 
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so that the effective potential is 

<!■(!/, () = ^ + (I) 

(10, 

The effective isothermal work performed in this reformu¬ 
lation is 


= ( 11 ) 

while the effective free energy change is Afeff = — 

with = $(?/,r) and = ^(y, 0). The Jarzyn- 
ski equality is 


correspond to the limit g —>■ 1, or kt —> 0 in this case. 
The work relations that we derive will refer to processes 
initiated when the system is in such a stationary state. 

The effective potential $(y,t) = ^niy) can be found 
in closed form for this model. We have 


‘ (l + lgr) 


- Kj' 


ksTo 

Kt 


In 1-1- 


KtX 

‘^-ksToJ 


V 


(17) 


From this we can write 




dtk ^ ° In 
Kt 



2kBTj ’ 


(18) 


exp (-VFeff[y]/A:B7o) = exp (-AF;ff/fcs?o), 

( 12 ) 

and all that remains is to cast the left hand side as an ex¬ 
pectation j dxP\x\ exp(—J^[kF]) in terms of a functional 
F of the rate of performance of work W{x,t). 


III. TSALLIS OSCILLATOR 


We illustrate how Eq. (12) can be interpreted for a 
model of 1-d particle motion in a harmonic confining po¬ 
tential and a temperature profile that varies quadratically 
with position measured from the centre of the potential. 
We consider (j){x,t) =(/>„ = K{t)x^f2 and 


Tr{x) = To 



ktx'^ \ 
^keTo) ’ 


(13) 


where kt is a parameter with the same dimensions as 
the spring constant k. This system is of interest since the 
Fokker-Planck equation for the pdf P{x, t) corresponding 
to the dynamics of Eq. (Q is 


so that the functional that appears in the Jarzynski 
equality is 


F[W] = 



dtW^^ 

KtX^ 


In 



ktx'^ \ 
2kBTo) ’ 


(19) 


using W = kx^ 12. This becomes clearer to interpret 
when written in the form 


The factor in square brackets is unity for an isother¬ 
mal bath, which can be demonstrated by writing = 
To -I- c{x), such that for small e we have [ro/(rr — 
To)] ln(T,/To) « (To/e)[e/To + 0[e^)] = 1 + 0(e) after 
which we take e —>■ 0, and the factor is less than or equal 
to unity for nonisothermal conditions. For an isothermal 
bath we would write F = W/kBTo- 

For a given x{t), the functional F of the associated 
work rate can be determined and (exp(—T)) set equal to 
exp(-[T($^(^)) - T($^(o))]/fcBTo) with 


In 


(Tr-To) 


Tr 


( 20 ) 


dP _ 9 f p ^ kB d{TrP) \ 
dt dx \m 7 TO 7 dx J ’ 

for which the stationary solution is 

® F (1 -b k/kt) a ktx'^ \~W~^ 

r(l-|-K/KT) V SfcsTo/ 

(15) 

namely a Tsallis distribution, proportional to a q- 
exponential function defined by 

e,{z) = (1 + (1 - , (16) 

here with q = [k + 2kt)/{k + kt) and 2 = —(k -|- 
kt)x‘^ /2kBTo. Such distributions have received atten¬ 
tion because of the ubiquity of experimental data that 
take this form |33| . They contrast with the more usual 
gaussians characteristic of equilibrium situations, which 


p:t{x) = 


Kt 


2TTkBTo 


exp[-T($«,)/fcBTo] cx 


= / dy 


dyexp(-$«,(y)/fcBTo) 
- —2k,Ikt — 1 


cosh[(KT/2fcBTo)^/^y] 

i r{l+ k/kt) 
^ T (1 -b k/kt) ’ 


( 21 ) 


noting that x = {2kBTol ktY^'^ smh[{KT/^kBToY^'^y] 
so that $K(y)/fcB'7o = 

[2k/kt + 1] lncosh[(KT/2fcBTo)^/^y]. In short, the 
Jarzynski equality for this system is 



To 

{Tr - To) 


In 



T (1 -b k{t)/kt) F (1 -b k{0)/kt) 
r (1 + i^{t)Ikt) r (1 -b k{0)/kt) ' 


( 22 ) 
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Figure 1. Illustration of a particle trajectory x{t) in a non- 
isothermal environment (indicated by background colour) as 
the spring constant changes in a sequence of steps. The ini¬ 
tial distribution of positions takes the stationary Tsallis form 
. At time r the distribution need not take the stationary 
form shown for comparison. 


Ki/ko 

Kt/ko 

g 

LHS of Eq. (25l 

1.01 

1.005 

1.5 

1.0000 ±0.0001 

1.05 

1.025 

1.49 

1.000 ±0.0011 

1.10 

1.05 

1.49 

1.000 ±0.0017 

1.20 

1.10 

1.48 

1.000 ±0.0021 

1.40 

1.20 

1.46 

1.001 ±0.0072 

1.75 

1.375 

1.44 

1.00 ±0.012 

2.0 

1.0 

1.33 

1.00 ±0.025 

4.0 

2.5 

1.38 

1.00 ±0.085 

6.0 

3.5 

1.37 

1.0 ±0.12 


Table I. Demonstration of the Jarzynski equality for the Tsal¬ 
lis oscillator for a range of work processes consisting of a step 
up and down in spring constant (from kq to ki and back 
again) separated by a period t = 4, with a variety of spa¬ 
tial temperature profiles specified by kt, in each case start¬ 
ing from a Tsallis distribution of position characterised by q, 
with m = y = ks = To = Ko = 1. Results are based on 
1000 numerical realisations for each set of parameters, using 
a timestep of 10“®. Errors in the mean are obtained assuming 
each realisation is an independent sample. The left hand side 
(LHS) of Eq. (|25[) is expected to equal unity. 



IV. JARZYNSKI EQUALITY FOR MULTISTEP 
PROCESSES 


where qi = 1 + Kr/iKi — Ki_i). This leads to the result 


We focus our attention on work processes consisting 
of a set of N abrupt changes in the spring constant, 
i.e. shifts Ki_i —> Ki at times t = ti for f = 1 to TV, 
with Kg = ^(0) ktv = The protocol of spring 

constant shifts is illustrated in Figure together with a 
representation of a stochastic trajectory x{t) in the pres¬ 
ence of a background temperature profile Tr{x). Tsallis 
distributions associated with stationary states at the 
beginning and end of the period are also shown: the sys¬ 
tem begins in a stationary state, but it need not end the 
process in such a condition. Note that a form of ergodic 
consistency holds here in the sense that the probability 
density is nowhere zero during the evolution: the stochas¬ 
tic dynamics can take the system from any initial position 
to any final position. 

For such a multistep process the functional P becomes 


N 


P{{Ki},{Xi}) = ^ 


, KtXI^ 


2=1 


K.'j' 


In 1 - 1 - 


2kBToJ ’ 


where Xi = x{ti). Hence 


(23) 


N 


exp(-j')= n ( 1+ 

2 = 1 
N 

n- 


2 \ -(Ki-Ki_l)/KT 


ktX^ 


(Kj - K^-i)a 


(24) 


/ N 

n< 


\ 2=1 


Am 

ksTo 


F (I -b KAf/Kr) F (1 -b Ko/kt) 


T (1 + km/kt) T (i -b kq/kt) ’ 

(25) 

where AWi = {Ki — Ki-i)xf /2 is the work associated with 
the ith step change in spring constant, such that the total 
work done is W — AWi- The notation now makes 
very apparent the resemblance to the isothermal Jarzyn¬ 
ski equality in Eq. Q. The brackets here refer to av¬ 
eraging according to (•) = /no^ {■)dXiPT{xo,Xi, ■ ■ ■ ,xn) 
where xq = a:(0) and pT is the probability that the dy¬ 
namics generate a discrete trajectory {xi}. If kt 0 the 
q-exponentials tend towards ordinary exponentials and 
the F functions can be represented using Stirling’s ap¬ 
proximation such that the right hand side of Eq. (251 
reduces to {kq/k^Y^^, as expected. 


It is straightforward to verify the identity (251 in sim¬ 


ple cases such as a process consisting of a step up from 
Kq to Ki at time = 0 followed by a step back down 
to K2 = Kq at t2 = T. Wc sclect initial positions ac¬ 
cording to the stationary distribution Pg^°{xo) given in 
Eq. (151 with a given kt, and solve the stochastic dy¬ 
namics in Eq. ([^ numerically in order to evaluate Xi 
and X 2 and hence AIFi and AW 2 - The expectation of 
e^j (—AkFi/fcsTo) Cqj (—Am/fc^Tb) over a number of 
realisations can be shown to be consistent with unity, 
within statistical uncertainty, as indicated in Table |T] An 
illustration of the distribution of n;ii e,. (-Am/fesTg) 
for a more complex cyclic process is shown in Figure 
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Fi|ure 2. Distribution of values of y = 

rii-i cyclic process with 

Ko = 1, Ki = 1.5, K 2 = 2, K 3 = 0.5, and K 4 = 1, 


with t 2 — ti = tz — t 2 = — tz = 0.1, and 

m = j = ks = To = kt = 1, obtained from 10® re¬ 
alisations with a timestep of lO”'*. The mean of the 
distribution is 0.999 with statistical uncertainty of 0.001, 
verifying the analogue Jarzynski equality (251 for this 
process. 



Figure 3. Probability distributions of work for process F con¬ 
sisting of a step up from kq = kt to ki = 2 kt (positive 
values of W), and process R specified by a step down from ki 
to Ko (negative values of W), both starting from a stationary 
state. The mean work for the F process is ksTo and for the 
R process it is —/cbTo/S. 


V. CROOKS RELATION FOR A STEP 
PROCESS 


Our next task is to determine the analogue of the 
Crooks relation for a Tsallis oscillator. The simplest 
demonstration involves an F process consisting of a step 
up in spring constant from kq to Ki > Kq starting from a 
stationary state, such that the corresponding R process 
is the step down from ki to kq. It is straightforward to 
evaluate the statistics of work performed from the statis¬ 
tics of position prior to the step change. We write 

pAW) = P:,°{xo)^, (26) 


with W = (ki — Ko)xQj2 > 0 such that dx^/dW = 
[2(ki — and so with appropriate normalisa¬ 

tion we have 


Pf[W) = 


Kt 


X e. 


r (1 -I- kq/kt) 
irksTo {ki - Ko)W\ T[^ + ko/kt) 

{kq + Kt)W 


QOF 


(ki - Ko)kBTo 


0(W), 


(27) 


with qqp = (kq + 2kt)/{ko + kt) and where 0 is the 
Heaviside function. Similarly we can determine the dis¬ 
tribution of work W = (kq ~ I‘^i)xq/2 < 0 for the reverse 
process starting from the stationary state with k = Kp. 


Pr[W) 


X e 


1 

tt/cbTo (ko - Ki)W _ 
f (ki -I- kt)W 


{Ko-ni)kBTo 


r (1 -I- kiIkt) 
r (i -b ki/kt) 

] e(-w), 


(28) 


with ( 7 oij = (ki-|-2k7’)/(ki-|-kt)- These two distributions 
are illustrated for the case of kq/kt = 1 and KijniT = 2 


in Figure The analogue Crooks relation is then 


Pf{W) 

PRi-W) 


r (I + Ki/K t) r (1 -|- Ko/ Kt) 

r (1 -|- Ki/ Kt) r (I -|- Ko/ Kt) 

X (l V^~ 

^ V ('^I “ Ko)A: bTo ) 

x(l- \ 

V (ko - ki)/cbTo/ 


(29) 


defined in this case for IF > 0, which may be written in 
the form 


Pf{W) 

Pr{-W) 


= exp(-[F($,J - F^($,„)]/fcBTo) 


X 



KjtW 

(ki - Ko)kBTo 


(ki-ko)/kt 


(30) 


and the final factor can be written as eqj^iW/kBTo) with 
qB = l- kt/(ki - Ko) or [eqj^{-W/kBTo)]~^ with qp = 
1 + kt/{ki — Kq). Once again, the notation makes it 
quite apparent that the Crooks relation for the pair of 
processes under isothermal conditions, in the form of Eq. 
is recovered when qp^R —)■ 1 and as kt —>■ 0. 

We can further conclude from these distributions that 
the mean work performed in the F process {W) p = 
f Wpp{W)dW is given by (ki — Ko)fcBTo/(2Ko — kt), as 
long as Ko > K7’/2, and the mean work for the R process 
is —(ki — Ko)/i:BTo/(2Ki — Kb), assuming that Ki > kb/ 2. 
For a process consisting of a step up, the establishment 
of a stationary state, and a step down, the mean work 
would therefore be 


{W)p + {W)r = 


2(ki - Ko)'^kBTo 
(2ko — kt)(2ki — Kt) ’ 


(31) 
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Figure 4. Verification of the analogue Crooks relation Eq. 
( |30| l for a step up from kq = kt to ki = 2kt and a step 
down from ki to kq, based on 10® particle trajectories with a 
timestep of 10“^, with m = 7 = fcs = To = kt = 1- 


which is never negative, and furthermore tends towards 
a known result |20] for such a process in the isothermal 
limit as kt 0. 

As an illustration of a possible experimental verifi¬ 
cation of the analogue Crooks relation, we have gen¬ 
erated trajectories according to Eq. Q for a step up 
from KqIkt = 1 to KiIkt = 2, and the corresponding 
step down, both starting from the appropriate station¬ 
ary state, in order to compute the work distributions. A 
plot of \nq^[pF{W/kBTo)/pR{-W/kBTo) exp{[F{^^^) - 
E($Ko)]/fcBTo)] against lE/fesTb should be a straight line 
with a gradient of unity, where In^ (z) = — l)/{l — q) 

is the g-logarithm |dd| . Such an outcome is apparent in 
Figure 


VI. PRACTICAL CONSIDERATIONS 


The work relations will depart significantly from the 
form taken under isothermal conditions when kt/h 1 
and we should consider whether current methods for es¬ 
tablishing a thermal profile in a particle trap can pro¬ 
vide such circumstances. Spring constants for particle 
trapping with optical tweezers in a liquid medium are 
typically of the order of 10“^ pN/nm HHun] and if we 
adopt this value for kt in Eq- (131, the thermal gradi¬ 
ent KRx/kB at a typical displacement of 10 nm from the 
trap centre would be about 10 K/nm. This is not an 
unattainable gradient in a solid system ciioni but could 
pose experimental challenges for a fluid system. 

In contrast, an approximately quadratic temperature 
profile over a spatial scale of microns with a peak ther¬ 
mal gradient of about 10 K//im can be established by 
laser illumination of a patterned nanostructure immersed 
in a fluid medium |29| . A particle can be confined 
within such a profile by thermophoresis (thermal diffu¬ 


sion) [15] or potentially by optical methods as well. Ther- 
mophoretic confinement provides much weaker trapping, 
with a spring constant of order 10“^ pN/nm. The effect 
is equivalent to the presence of a potential given in di¬ 
mensionless form by (j){x)/kBT = ST^Tr{x), where T is 
an effective constant ambient temperature and St is the 
Soret coefficient that compares the strength of thermal 
diffusion, in response to a temperature gradient, with 
that of ordinary diffusion down a density gradient. For 
a quadratic temperature profile the confinement is har¬ 
monic with an effective spring constant n = TStKt and 
strongly nonisothermal conditions will therefore be es¬ 
tablished if TSt ~ 1- Such a thermophoretic trap was 
demonstrated in |2^ using 200 nm polystyrene spheres 
with St ~ 3 K“^, and hence kt 10“^k, but by trap¬ 
ping smaller species the Soret coefficient could potentially 
be reduced to 10“^ K~^. There remains the challenge 
of ensuring that such a trap is strong enough to retain 
a particle for a long enough time to perform measure¬ 
ments, and compromises will need to be made between 
the desired temperature profile and such requirements. 
Nevertheless, it seems possible that with further techno¬ 
logical development the conditions of interest might be 
attainable. 


VII. CONCLUSIONS 

Work relations exist for systems that are maintained 
away from thermal equilibrium by constraints, and a 
number of such cases have been discussed |151 - H5] . The 
situation we have explored here is of interest since the 
system in question, when described using overdamped 
dynamics, is governed by Tsallis statistics in its station¬ 
ary state. Such statistical properties have been noted 
in a variety of physical systems |33|, and contrast with 
the more usual gaussian statistics of thermal equilibrium. 
The analogues of the Jarzynski equality and Crooks rela¬ 
tion for such a system, constructed for certain processes, 
take a particularly transparent form when we make use of 
the g-exponential functions associated with Tsallis statis¬ 
tics. 

Stochastic entropy production can provide the con¬ 
ceptual underpinning of work relations in nonisothermal 
as well as isothermal conditions, employing the usual 
Gibbs-Shannon form rather than the Tsallis entropy. We 
believe that experimental verification of these identities 
might be possible in small systems where the appropri¬ 
ate nonisothermal environment can be established and 
controlled. 
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